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SUMMARY 


This  thesis  deals  with  the  investigation  of  a  sampled- 
data  system,  with  a  quantizer  in  the  forward  loop. 

Two  methods  of  simulating  the  nonlinear  element  (quan¬ 
tizer)  were  successfully  developed.  Both  utilized  analog 
computer  components  which  consisted  mainly  of  comparators 
and  amplifiers. 

The  transient  response  of  an  error-sampled,  second- 
order,  unity-feedback,  quantized  system  was  studied.  Treat¬ 
ing  the  quantizer  as  a  nonlinear  element,  the  state-variable 
technique  was  used  for  the  purpose  of  analysis.  Calculated 
results  (by  means  of  a  digital  computer)  and  results  ob¬ 
tained  on  the  analog  computer  were  in  agreement. 

A  stability  criterion  by  Tsypkin  ,  based  on  the 

work  of  Popov,  was  used  to  determine  the  globally  asymp¬ 
totic  stability  of  the  system  with  respect  to  the  gain  of 
the  linear  plant.  All  information  in  the  frequency  domain 
was  converted  into  a  gain-phase  plot,  from  which  the  stab¬ 
ility  criterion  could  easily  be  interpreted. 

A  recently  modified  criterion  by  Jury  and  Lee  was 
also  considered  and  it  yielded  much  less  conservative  re¬ 
sults.  In  this  case,  a  digital  computer  was  used  to  apply 
the  criterion. 

Another  study  carried  out  in  the  Z-domain  indicated 
that  a  filter  was  effective  in  compensating  the  system  for 
both  cases  so  that  a  larger  gain  was  possible;  or,  in  other 


aye 


words ,  the 
The  filter 


relative  stability  of  the  system  was  increased, 
so  found  was  physically  realizable. 
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1 .  Introduction 

In  sampled-data  control  systems,  the  digital  com¬ 
puter  is  becoming  more  and  more  important  because  of  its 
high  accuracy,  speed,  versatility  and  flexibility.  From 
the  inherent  nature  of  numerical  manipulation,  the  problem 
of  quantization  arises  naturally.  Especially,  in  dealing 
with  hybrid  simulation  or  computation,  in  the  areas  of 
digital  to  analog  conversion  or  vice  versa,  the  process 
requires  the  unavoidable  rounding  off  of  numbers  which  can 
be  considered  as  quantization  errors.  When  the  accuracy 
of  the  analog  device  is  specified,  the  capacity  of  the 
digital  device  in  the  sense  of  word-lengths  is  determined 
accordingly  in  terms  of  the  fineness  of  quantization. 

What  is  more  important  is  the  problem  of  stability,  when 
there  is  a  quantizer  in  the  system. 

This  thesis  deals  with  the  stability  of  a  sampled- 
data  system  with  a  single  quantizer  in  the  forward  loop. 

The  general  approach  is  to  regard  the  quantizer  as  a 
nonlinear,  time-invarient ,  memoryless  element,  which  is, 
in  a  sense,  a  relay  with  multiple  levels  of  output.  Thus, 
the  system  has  become  a  nonlinear  sampled-data  system 
which,  due  to  lack  of  completeness  and  uniqueness  of  basic 
principles  in  this  area,  causes  much  difficulty  in  analysis 
and  synthesis.  Moreover,  as  a  result  of  the  difficulty  of 
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describing  the  quantizer  in  mathematical  terms,  the  gen¬ 
eral  approach  to  the  analysis  of  this  kind  of  situation 
is  by  graphical  methods «  There  still  remains  the  problem 
of  how  to  simulate  a  quantizer  and  compensate  for  stab¬ 
ility  of  such  systems.  In  this  thesis,  a  new  method  for 
solution  of  these  problems  is  presented. 
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2 .  Analog  Simulation  of  a  Quantizer 

A  quantizer  is  a  nonlinear,  time-invariant,  memory¬ 
less  element  with  an  input-output  transfer  characteristic 


input 


At  the  present,  there  are  several  methods  available 
for  constructing  a  quantizer  > (3) a  However,  none 


of  these  is  satisfactory  because  of  the  following 
disadvantages;  excessively  long  response-time,  complex 
solid-state  circuitry,  inaccuracy,  drift  of  output  as 
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time  elapses,  wear  and  tear  of  mechanical  parts,  and  very 
limited  number  of  steps  achieved. 

Using  the  analog  computer,  a  very  simple  technique 
was  found  to  be  satisfactory.  It  is  referred  to  as  the 
"Parallel  Method"  of  simulating  a  quantizer.  This  is 
implied  by  the  parallel  arrangement  of  the  components  used. 
On  the  other  hand,  as  a  conjugate  case,  a  "Series  Method" 
was  also  found  possible. 

2(A)  Parallel  method  for  simulating  the  quantizer 

A  block  diagram  illustrating  the  parallel  construc¬ 
tion  of  the  quantizer  is  given  in  figure  2. 


Figure  2.  Block  diagram  showing  parallel  method  of 
simulating  a  quantizer. 


■ 
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Comparators,  A,  B,  C,  D  „  .  .  ,  are  illustrated  in 
figure  3 . 


Figure  3.  Block  diagram  of  a  double-pole- 
double-throw  comparator . 

Inputs  and  outputs  are  indicated  by  direction  of  arrows. 
The  function  of  the  comparator  in  figure  3  is  illustrated 
by  table  1. 


Algebraic  sum  of 

inputs 

el  +  e2 

outputs 

e01 

e02 

>0 

e3 

e4 

<  0 

e5 

e6 

Table  1.  Truth  table  of  the  comparator. 

The  dotted-line  portion  of  the  block  diagram  shown  in 
figure  2  gives  the  first  quadrant  of  the  transfer  curve, 
and  another  comparator  is  used  to  change  sign  when  the 
input  is  negative. 


' 
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Operation  of  the  scheme  in  figure  2  can  be  understood 

better  by  considering  two  typical  examples  as  follows.  If 

3L  jL  13 

the  input  is  between  JI  and  JI  ,  output  of  comparator 

D  is  4q,  and  that  of  C  is  2q.  The  switching  value  of  B 
1  13 

is  now  -  jl+4J[  +  2^  =  —  Ji.  The  net  quantizer  output 

3b  3  1 5 

is  4q  +  2q  =  6q.  If  the  input  is  between  -y  1  and  — —  X , 
outputs  of  D  and  C  are  still  respectively  4q  and  2q.  How¬ 
ever,  input  of  the  first  one,  B,  can  now  overcome  its 

13  ft 

switching  value  of  —  JI ,  thus  giving  lq  as  its  output. 

The  net  quantizer  output  is  now  6q  +  lq  =  7q.  The  first 
comparator,  A,  is  used  to  yield  the  same  sign  for  both 
the  input  and  output.  Following  the  same  steps  as  above, 
the  number  of  steps  possible  and  the  number  of  comparators 
used  is  tabulated  as  follows: 


*  Number  of 

Comparators  Used 

Total  Number  of 
Steps  Possible 

1 

2 

2 

6 

3 

14 

4 

30 

5 

62 

6 

126 

7 

254 

8 

510 

9 

1022 

10 

2046 

11 

4094 

12 

8190 

13 

« 

16382 

• 

• 

* 

• 

• 

Table  2.  Relation  between  number  of  comparators  and  number  of 
steps . 

*  For  each  case,  one  more  comparator  has  to  be  used  for  sign- 
inversion  . 


«■  . :  r  so  S 

30  S'  1  'O 

.  .4  -  !  .  ’  ■ . 

i .. 

• 

.  i  - 

*  I-  -  >  •$  - 

'*•  .  .  i  .• ' 

e  .  ■  '  ■  '  •;  ri'4 


•  ■  ■;  r.  '  . 

£ 

' 

L  <  ' 

,  .  ■  ■-  -  :  •  '  ■  ' 
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A  schematic  diagram  of  254  steps  (+  127  steps) ,  each 
of  which  is  one  volt,  is  shown  in  figure  4.  For  the  pur¬ 
pose  of  convenience,  the  same  scale  (q  =  Jl  =  1)  is  used 
for  both  the  input  and  output.  This  scheme  requires  eight 
comparators,  fifteen  amplifiers  and  fourteen  pots.  For 
a  computer  with  one  hundred  or  more  amplifiers,  ten  or 
more  comparators,  one  hundred  or  more  pots,  this  does  not 
create  any  difficulty  due  to  a  shortage  of  components. 


Figure  4. 


Schematic  diagram  of  parallel  method  (q  ~  l  -  1) 
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2(B)  Series  method  of  simulating  the  quantizer 


A  series  method,  a  conjugate  of  the  parallel  type, 
is  also  possible.  This  is  shown  in  figure  5. 


Figure  5.  Block  diagram  showing  series  method 
for  simulating  a  quantizer. 

9 

An  extra  comparator,  A  ,  is  also  necessary  for  proper 
consideration  of  signs  of  input  and  output. 

The  quantizing  process  is  again  illustrated  by  two 
typical  examples.  If  the  input  is  between  &  and  •—  A, 
comparator  B  yields  4q.  Input  to  C  is  ein  -  4  JL  ^  X 
which  causes  its  output  to  be  2q.  Input  to  D  is  e^n  -  4  X 
2  Ji  4  which  is  not  sufficient  to  cause  any  output.  The 
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net  result  is  6q.  If  the  input  is  now  between  ™-  X  and 
15  '  « 

— 2*  X,  outputs  from  B  and  C  are  respectively  4q  and  2q. 
Input  to  D  is  e^n  -  4  J  -  2  j(  ^  which  causes  an  output 
of  lq„  The  net  result  is  7q„  The  relation  between  the 
number  of  comparators  necessary  and  the  steps  feasible  was 
found  to  be  exactly  the  same  as  that  given  in  table  2 
of  the  previous  case. 

A  complete  schematic  diagram  with  +  127  steps  of 
unity  scale  is  shown  in  fig.  6.  This  scheme  yields  id¬ 
entical  results  as  the  parallel  type.  The  total  number 
of  computer  components  used  is  equal  to  that  of  the  previous 


case . 
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Figure 


5.  Schematic  diagram  of  series  method  (JL  =  q  =  1) 
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2(C)  Merits  of  both  methods 

As  in  the  case  of  the  ideal  sampler  or  relay,  the 
construction  of  an  ideal  quantizer  is  a  physical  impossib¬ 
ility.  The  merits  of  the  quantizer  realized  by  the  above 
two  methods  depend,  of  course,  to  a  very  great  extend  on 
the  components  used;  i.e.,  the  comparators,  amplifiers, 
and  the  accuracy  by  which  the  pots  are  set.  The  compara¬ 
tors  available  are  electromechanical  devices  with  a  dead 
zone  of  +  10  m.v.  and  the  response-time  of  two  milliseconds 
maximum.  Pot-setting  accuracy  as  indicated  by  the  digital 
voltmeter  on  the  computer  is  also  limited  to  the  order  of 
10  m.v.  The  actual  accuracy  of  the  quantizer  is  satis¬ 
factory  for  normal  usage  on  the  computer.  However,  when 
the  maximum  amplitude  of  the  overall  response  of  the 
system  does  not  exceed  one  hundred  volts  as  limited  by 
the  rated  voltage  of  the  computer,  the  quantizing  steps 
can  be  scaled  up  for  even  greater  accuracy.  For  example, 
if  the  maximum  amplitude  of  the  response  of  the  system  is 
not  greater  than  twenty  volts,  a  one-volt-step  quantizer 
can  be  scaled  up  to  one  having  five- volt  steps;  which  means 
that  the  other  amplitudes  have  to  be  correspondingly  scaled 
up.  In  this  case,  accuracy  is  increased  by  a  factor  of  five. 

Experiments  performed  indicated  that  any  correctly 
quantized  level  in  the  output  could  be  reached  in  less 
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than  approximately  five  milliseconds.  For  a  sampled-data 
system  with  a  sampling  period  of  the  order  of  one  second 
or  even  one-tenth  of  a  second,  the  duration  of  five  milli¬ 
seconds  is  negligible  in  most  experimental  work. 

A  great  advantage  of  these  two  methods  over  all  others 
is  that  the  desirable  amount  of  each  step  can  be  adjusted 
very  easily;  i.e.,  the  values  of  q  and  ji  are  changed 
merely  by  individual  pot-setting.  In  addition,  there  is 
no  drift  of  voltage  in  the  output  as  normally  occurs  in 
analog  memory  devices.  Therefore,  the  methods  discussed 
are  satisfactory  in  terms  of  flexibility,  large  number 
of  steps  possible,  fast  response-time,  convenience,  sim¬ 
plicity  and  cost  of  construction. 


bnoo^e  £«no  v&b-to 
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3  Transient  Response  of  S.D.  System  with  a 
Quantizer  in  the  Forward  Loop 


A  second-order ,  error- sampled ,  unity-feedback  system 
was  chosen  for  transient  response  analysis  and  it  is  shown 
in  figure  7. 


The  unit  for  the  quantizer  is  unity  (-Jt™  q  ®  1)  on  both 


puter,  experiments  were  performed  with  various  initial 
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conditions  and  gains  of  the  linear  plant .  Phase-plane 
diagrams  were  used  to  observe  the  transient  response. 

With  different  values  of  gains,  stable  response  or  limit- 
cycle  oscillation  was  observed  (the  critical  value  of 
gain  for  stability  will  be  discussed  later) . 

The  state-variable  technique  was  adopted  for  theor¬ 
etical  consideration.  The  quantizer  is  treated  as  a  variable- 
gain  amplifier  which  has  a  constant  gain,  K(KT) 
during  the  interval  of  each  sampling  period.  The  signal 
flow  graph  is  shown  in  figure  8„ 


From  the  signal  flow  graph  in  figure  8,  where  =  C, 
x2  “  *1  ' 


x1  (s) 


(s+1) - 


x^(KT)  + 


1 

s  ( s+1) 


x2(KT)  + 


-  KKk  ~ 
-s^  ( s+1)  - 


r  (KT) 
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r  -KK. 


x2(s)  = 


k 


s  ( s+l) 


x  ^ ( KT )  + 


s+l 


x2(KT)  -f 


KK. 


k 


s  ( s+ 1 ) 


r  (KT) 


taking  the  inver se-Laplace  transform; 


x. 


(k+1)  T 


1  -  KK,  (T  -  1  +  e  T) 
k 


xx(KT)  +  (1  -  e  T)  x 


+  KkK  (T  -  1  +  e“T)  r (KT) 


x. 


(K+l)  T 


-K^K  (1  -  e"T)  xx(KT)  +  e”T  x2 (KT) 


— KT 

+  K^K  (1  -  e  )  r  (KT) 


in  the  form  of  matrix  algebra; 


X 


(K+1)T 


=  A  (T)  X (KT)  4-  B  (T)  R  (KT) 


for  intersampling  response,  for  0  <  A  <  1 


X 


(K+A) T 


=  A (AT)  X (KT)  +  B (AT)  R  (KT) 


In  our  case,  let  T  =  1  sec.,  r(t)  =  0,  the  state- 
transition  equations  become; 


x 


1 


(K+A) T 


L 


-AT 

1  -  K^.K  (AT  -  1  +  e  ) 


x1(KT)  +  (1  -  e  AT) 


x  2 ( KT )  ... 


2(kt) 


(1) 
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x 


2 


(K+A) T 


-K^Kd  -  e"AT)  xx(KT)  +  e  AT  x2  (KT) 


(2) 


For  manipulation  of  the  transient  response  point  by 
point,  a  digital  program  was  written  in  FORTRAN  4  and  the 
solution  obtained  on  the  IBM  7040 „  The  complete  program 
is  given  in  Appendix  1.  The  quantizer  in  the  digital  pro¬ 
gramming  was  formulated  by  use  of  the  relation  between 
floating-point  and  fixed-point  numbers.  The  result  so 
obtained  checked  well  with  that  obtained  on  the  analog 
computer . 

The  usefulness  of  the  state-variable  technique  with 
the  aid  of  a  digital  computer  is  certainly  not  limited  to 
a  second-order  system.  With  systems  of  higher  order  and 
nonunity  feedback,  the  transient  response  can  still  be 
computed  step  by  step  in  the  same  fashion. 
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4  Stability  Criterion  (Tsypkin's  Case) 

Based  on  the  work  by  Popov  for  the  case  of  a  nonlinear 

( 11) 

continuous  system,  Tsypkin  has  derived  a  criterion  for 

the  determination  of  absolute  stability  in-the-large  for 
nonlinear  sampled-data  systems „  The  basic  system  is  shown 
in  figure  9. 


Figure  9*  Basic  system  for  Tsypkin's  case. 

If  0(0)  =0,  0  <C  <1  k  for  of  0,  and  if  the  linear 

plant  F(Z)  is  stable;  then  the  autonomous  system  is  glob¬ 
ally  asymptotically  stable  provided  that  the  following 
condition  is  satisfied? 

™  +  Re  F(Z)  —  0  for  every  value  of  Z  on  the 
unit  circle  jz|  =1  ;  where  k  is  the  sector  to  which  the 
nonlinear  characteristic  is  confined  as  shown  in  figure  10. 


srfi  no  S  eulsv  yiavs  toi  0 
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Figure  10.  The  nonlinear  sector  k. 

The  condition  that  requires  0  <  k  for  a  f  0 

can  be  modified  if  a  quantizer  is  considered  to  be  the  non¬ 
linear  element,  since  it  does  not  hold  for  the  range  where 
J a |  This  can  be  done  by  allowing  the  steady-state 

error  to  be  within  the  range  of  +  -^instead  of  zero^^  . 

The  result  is  the  same  since  the  output  of  the  quantizer 
is  zero  within  this  range. 

The  question  of  absolute  stability  has  thus  become 
considerably  easier  as  it  refers  only  to  the  frequency 
response  of  the  linear  plant  against  the  region  of  con¬ 
finement  of  the  nonlinearity;  and  this  is,  as  the  word 
"absolute"  implies,  independent  of  the  shape  of  the  trans¬ 
fer  curve  of  the  nonlinearity.  However,  this  is  a  suf¬ 
ficient  but  not  necessary  condition,  since  the  result  so 


«?:  sr  s  non  t. 
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obtained  is  quite  conservative.  This  criterion  will  be 
investigated  before  the  modified  form  is  considered. 

A  study  was  carried  out  to  determine  the  maximum 
gain  of  the  linear  plant  for  stability  once  the  sector 
of  confinement  of  the  nonlinearity  is  fixed  by  the  specific 
element  chosen.  In  the  quantizing  case  studied,  the  sector, 


The  idea  of  frequency  response  implies  the  commonly 
used  polar  diagram.  Since  the  maximum  gain  of  the  plant 
is  required,  a  minimum  value  of  the  real  part  of  the 
frequency  response  is  expected  when  the  gain  is  set  at 
unity;  or 


min  Re 


F(Z) 


(4) 


for  every  point  on  |  Z  I  =  1. 

The  above  condition  yields  the  polar  diagram  of  a  second- 
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order  system  shown  in  figure  12. 


Figure  12.  Polar  plot  for  Tsypkin ' s  criterion 
of  a  second-order  system. 

Since  the  frequency  response  is  obtained  by  having 
Z  take  on  values  on  the  unit  circle  in  the  Z-domain;  and 
the  frequency,  co,  appears  as  an  exponent  in  every  term; 
there  is  no  quick  and  systematic  way  to  obtain  the  response 
as  could  be  done  with  a  Bode  plot  in  the  continuous  case. 
The  tedious  job  could  be  handled  very  easily  by  a  digital 
computer,  and  the  minimum  real  part  obtained  directly. 

A  new  graphical  method  is  also  possible. 

A  commonly  used  system,  error-sampled  and  followed 
by  a  zero-order-hold  before  the  nonlinearity  is  certainly 
not  of  a  form  to  which  Tsypkin ' s  criterion  can  be  directly 
applied.  The  difficulty  could  be  solved  by  interchanging 
positions  of  the  appropriate  elements  in  the  system  as 


.. 


■ 
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' 
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shown  in  figure  13. 


Figure  13.  A  new  location  for  the  quantizer. 

It  can  be  shown  that  the  two  schemes  are  identical  as  far 
as  input  to  G(s)  is  concerned,  so  that  Tsypkin's  criterion 
can  be  directly  applied. 

The  linear  plant  is  now  F(s) 


K 

s  ( s+1) 


and 


F(Z)  =  ^ 


F  ( s) 


(Z-l) (Z-e"1) 


for  T 


second,  K  =  1, 


where  p  = 


-  (1  -  2e  x) 


0.717 


For  this  case,  results  obtained  on  the  digital  computer 


indicated  that  the  minimum  real  part  was  “1.5 


tical  result  is 


<!+ 1) 


the  analy- 


.  For  k  =  2,  and  substituting 


the  numerical  values  into  equation  4  results  in 


1 

2 


4  -1.5K 


or 


K  4 


and 
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the  maximum  gain  is  found  to  be  1/3. 

In  addition  to  the  above  method,  a  graphical  tech¬ 
nique  is  also  possible  for  interpreting  the  criterion  in 
a  different  manner  which  enables  the  problem  of  compen¬ 
sation  to  be  handled. 

All  information  in  the  frequency  response  is  now 
converted  into  a  gain-phase  plot  with  the  amplitude  in 
decibels  for  convenience,  and  the  phase  angle  in  degrees. 
The  line  which  becomes  a  curve  in  the  gain-phase  plot 
will  be  referred  to  as  the  "Tsypkin  curve"  and  is  shown 
in  figure  14.  In  this  case,  for  the  condition  to  be 
satisfied,  the  Tsypkin  curve  has  to  be  above  or  just  tan¬ 
gent  to  the  curve  plotted  for  the  linear  plant  for  all 
frequencies.  Therefore,  the  Tsypkin  curve  is  plotted  on 
transparent  paper;  and,  in  order  to  maintain  the  same 
phase  relation,  it  is  shifted  up  and  down  on  the  gain- 
phase  plot  until  the  Tsypkin  curve  is  tangent  to  that  of 
the  linear  plant.  The  intersection  of  the  Tsypkin  curve 
with  the  vertical  axis  indicates  the  minimum  value  of  the 
real  part  of  the  frequency  response  of  F(Z). 

For  example,  the  second-order  system  shown  in  figure 
13  is  considered  again.  The  frequency  response  of  such 
a  system  is  plotted  in  figure  15  in  which  the  Tsypkin 
curve  is  superimposed.  The  intersection,  point  A,  gives 
a  value  of  approximately  3.75  db  or  1.54  (actual  value 
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being  1.50).  The  error  is  about  2.6%  which  is  considered 
quite  satisfactory  for  a  graphical  method. 

In  conclusion,  this  approach  gives  a  new  method  of 
obtaining  the  required  stability  information  graphically 
with  reasonable  accuracy  as  far  as  absolute  stability  is 
concerned.  Also,  it  proves  quite  useful  for  compensation 
as  will  be  discussed  in  section  6. 


criterion  of  a  second-order  system. 
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5  A  Phase- lead  Filter  in  the  Z-domain 


The  properties  of  a  phase-lead  filter  will  be  dis¬ 
cussed  and  the  conditions  sufficient  for  compensation  will 
be  derived  to  facilitate  compensation  in  the  Z-domain. 

To  begin,  a  very  simple  filter  with  a  real  pole  at 
B,  and  a  real  zero  at  A,  as  shown  in  figure  16,  is  assumed 
in  the  Z-domain;  where  B  >  A  and  B  1,  A  <C  1 . 


Figure  16.  Pole-zero  configuration  foi  the  lead-network 


Hence  the  transfer  characteristic  is 


D(53)  = 


Z-B 

Z-A 


1-A 


Such  a  filter  has  a  minimum  amplitude  of  — -  at 


<oT 


0 


H'A 

and  a  maximum  value  of  "rip  at  oYT 
a  maximum  value  of  lead  of  d 


—  sin 

max 

_  *  -1  a+b 

at  the  frequency  of  T  -  cos  i’+Xb  * 


77 .  The  phase  reaches 
1  B-A 


I  AB  occurs 

These  relationship 


are  plotted  in  figure  17. 


' 


, 


' 
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Figure  17 0  Frequency  response  of  the  lead-network „ 
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6.  Compensation  (Tsypkir^s  Case) 


The  frequency  response  of  a  linear  second-order  plant 
is  plotted  in  figure  18  and  is  depicted  by  the  solid  line, 
GA\^  (db) 
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By  incorporating  the  properly  chosen  lead-network, 
the  original  gain-phase  plot  is  shifted  down  in  amplitude 
for  relatively  constant  phase.  This  allows  for  a  larger 
gain  of  the  linear  plant,  since  the  Tsypkin  curve  1  can 
be  shifted  until  it  is  tangent  to  the  new  gain-phase  plot; 
thus  giving  a  smaller  minimum  real  part  of  the  compensated 
system  (N  <  M) . 

All  work  is  carried  out  in  the  Z-domain  and  the  in¬ 
formation  is  obtained  from  a  gain-phase  plot  without  any 
further  transformation.  The  use  of  the  ordinary  Bode 
plot  in  this  case  is  found  to  be  impossible. 

For  the  convenience  of  design  of  a  proper  lead- 
filter,  the  following  relationships  are  necessary? 

1  +  sin  (<t>  -  cd  T) 

B  .  rn  m 

sin  4>  +  cos  cd  T 

m  m 

sin  -  cos  cd^T 

A  -1  +  sin  “  a)  T) 

m  m 

where  B  and  A  are  locations  of  the  pole  and  the  zero  on 

the  real  axis  inside  the  unit  circle. 

Graphs  showing  the  frequency  response  of  the  filter 

are  plotted  in  figures  19  to  26  for  various  values  of  ’a)  T 

and  <t>  which  determine  the  locations  of  pole  and  zero  of 
m 

the  filters  to  be  used  as  compensators. 
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In  the  following  autonomous  system,  shown  in  figure  27(a), 
Fc(Z)  includes  both  the  original  plant  and  the  lead-filter. 


Figure  27(a).  Compensated  system. 


Figure  27(a)  is  identical  to  figure  27(b)  for  analysis 
and  synthesis. 


Figure  27 (b) .  Compensated  system, 


In  figure  27(b),  D  (s)  is  shown  to  contain  a  zero-order-hold, 


The  linear  plant 


G ( s)  and  ZOH 


is  plotted  in  figure  28 


which  yields,  as  before,  a  minimum  real  part  of  -1.5  with 

1  1 
K  =  -j.  If  the  linear  plant  is  to  have  a  gain  of  it 

implies  that  the  Tsypkin  curve  should  be  shifted  down  to 
yield  a  gain  of  zero  db  at  -180°  phase  shift  (indicated  by 
the  dotted  line  in  the  same  figure) .  Considering  the  orig¬ 
inal  gain-phase  plot  of  figure  28,  the  curve  should  be  modi¬ 
fied  so  that  between  the  range  of  frequencies  o)^  =  0.664 
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and  a)^  =  0,  the  Tsypkin  curve  should  be  above  the  gain-phase 
plot  for  all  phase  angles  in  order  to  maintain  absolute 
stability.  Judging  from  figures  19  to  26,  by  trial-and- 
error,  the  following  optimal  values  are  picked? 


$  =  10° 
m 

cd  =0.2 

m 

which  yield  A  =  0.78639,  B  =  0.84469. 
Therefore,  the  lead-filter  is 


D(Z) 


Z  -  0.84469 

Z  -  0.78639 


The  compensated  system  is  also  plotted  in  figure  28. 

This  result  is  close  to  the  required  gain  of  the  plant. 

In  fact,  as  checked  by  the  digital  computer,  the  exact 
value  of  the  minimum  real  part  is  -0.9755.  The  compen¬ 
sated  system  is  now 

(  *  =  0 „ 368K (Z  +  0.717)  (Z  -  0.84469). 
cl  '  (Z-l)  (Z  -  0.368)  (Z  -  0.78639) 

where  K  =  ~  . 

At  this  point,  it  may  be  appropriate  to  clarify  the 
relative  locations  of  the  various  components  in  the  system. 
The  original  system  without  compensation  is  shown  in 
figure  29. 
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Figure  29.  System  without  compensation. 


To  convert  the  block  diagram  into  a  form  where  the  stabi¬ 
lity  criterion  can  be  directly  applied  without  changing 
the  system,  the  scheme  shown  in  figure  30  resulted. 


h 


<2> 


a 


ZoH 


GCs) 


f(s)- 


Figure  30.  A  form  where  the  stability  criterion 
can  be  directly  applied. 

The  new  linear  plant  is  one  that  contains  both  the  ZOH 
and  G(s).  On  the  basis  of  this  linear  part  the  critical 
gain  is  derived  for  absolute  stability.  No  change  is  made 
on  the  system.  For  compensation,  a  lead-network  is  de¬ 
signed  in  the  Z- domain .  In  order  to  maintain  the  form 
of  the  original  error-sampled  system  followed  by  a  ZOH, 
and  in  order  that  the  system  without  compensation  consists 
of  a  ZOH  and  G(s),  the  compensating  filter  D(Z)  should  con¬ 
tain  a  ZOH  in  front.  The  open-loop  system  is  given  in 
figure  31. 


' 
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Figure  31.  Open- loop  system  with  compensation. 

In  the  system  shown  in  figure  31,  the  arrangement  of  the 
block  diagram  is  made  for  convenience  in  theoretical  work 
and  the  actual  open-loop  system  takes  on  the  form  shown 
in  figure  32. 


Figure  32.  Actual  open-loop  system  after  compensation. 


X(Z)  is  the  effective  compensation.  Considering  the 

Z— B 

scheme  in  figure  31  where  D(Z)  =  ^  (B  >  A)  ,  it  is  found 

to  have  the  arrangement  shown  in  figure  33.  Note  that 
D (Z)  f  X(Z)  . 


ZoH 


p(3) 


V 

T 


r 


D(2> 


Figure  33.  Block  diagram  of  compensator. 


4 


43 


D(s)  can  be  synthesized  by  ordinary  methods  of  realization, 
since  the  number  of  poles  is  equal  to  the  number  of  zeros,  and 
the  pole  and  zero  are  both  simple,  real  and  inside  the  unit 
circle  in  the  Z-plane.  The  zero  is  always  positive  as  is 
the  pole  in  most  cases.  When  the  pole  is  positive,  it 
is  possible  to  realize  B{s)  by  a  series  pulsed-data  net¬ 
work  or  a  digital  computer.  Otherwise,  if  the  pole  is 
negative,  D(s)  may  be  realized  on  a  digital  computer. 

Considering  the  example  where 

D(Z)  =  ,  it  can  be  realized 

for  the  series  case  shown  in  figure  34. 


D  OO  = 

2.4 


DU) 


Figure  34.  Realization  by  a  series  pulsed-data  network. 

For  the  feedback  pulsed-data  network,  it  is  always  pos¬ 
sible  for  a  realizable  zero,  however,  in  the  actual  arrange¬ 
ment  of  the  system,  the  ZOH  cannot  be  separated  from  D(s) . 
Therefore,  this  method  is  disregarded. 

For  the  time-delay  digital  programming  method,  in 
order  to  separate  a  ZOH  from  D(Z) ?  it  can  be  shown  that 
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D(Z) 


Z-B 

Z-A 


D  (Z) 


D  (Z)  . 


This  scheme  is  shown  in  figure  35  where 
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Figure  35 0  Realization  by  a  digital  computer. 


Figure  35  is  redrawn  below  in  figure  36(b)  showing  the 
digital  program  inserted,.  Figure  36(a)  shows  a  unit  time- 
delay  which  is  part  of  the  scheme  in  figure  36(b) „ 


Figure  36(a).  A  unit  time-delay  of  T  seconds. 


Figure  36(b).  Realization  by  a  digital  computer. 
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Addition  of  a  ZOH  and  a  sampler  in  front  of  the  digital 
computer  does  not  change  the  input  to  G(s) ,  since  the  out™ 
put  of  the  digital  computer  must  also  be  sampled  and  held. 


*< 

7  Stability  and  Compensation 

(Jury  and  Lee  * s  Case) 
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A  modified  criterion  based  on  Tsypkin ' s  result  has 

(4) 

been  derived  by  Jury  and  Lee  .  They  have  found  a  suf¬ 
ficient  condition  of  absolute  stability  for  this  kind  of 
system  which  is  much  less  conservative. 

If,  for  the  transfer  curve  of  the  memoryless  non- 


the  nonlinear  gain  function  is  assumed  continuous,  and 
satisfies  the  conditions? 

<t>(0)  =0,  k  >  >0  for  o  f  0,  j^0)  )<  k'  , 

II 

or  k  (  k  where  k  is  the  local  slope  of  the  nonlinearity  ? 

the  linear  part  must  be  stable,  and  lim  G(s)  =0. 

»  s-«~  oo 


*  in  the  case  where  the  nonlinear  function  is  piece-wise- 
linear,  the  condition  for  the  local  slope  is  satisfied  if 
the  slope  of  each  linear  segment  is  less  than  or  equal  to  k 
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The  theorem  states  that;  such  a  nonlinear  sampled-data 
system  is  absolutely  stable  if  a  non-negative  q  exists 
such  that  the  following  condition, 


47 


Re  F (Z) 


1  +  q(Z-l) 


t 


2 


(Z-l)  F ( Z) 


2 

£  0 


is  satisfied  on  the  unit  circle  Z  =  e^T  .  .  .  (a)  .  Again 
the  same  assumption  as  before  is  made  that  the  condition 
k  >  — >  0  for  o  f  0  is  relaxed  for  finite  steady-state 
error  with  magnitude  less  than 

It  is  obvious  that  by  letting  q  =  0,  this  theorem 
is  reduced  to  the  previous  Tsypkin's  criterion  which  can 
therefore  be  considered  as  a  special  case  of  this  one. 
Examples  by  Jury  and  Lee  are  shown  in  some  instances  to 
yield  results  of  approximately  twice  the  gain  as  compared 

to  those  found  by  Tsypkin's  case  for  the  same  system. 

( 5) 

In  a  later  paper  by  the  same  authors  ,  the  above 
condition  has  been  extended  to  some  specific  subclasses 
of  the  basic  system.  One  of  the  extended  cases  which  is 
suitable  for  application  to  systems  containing  a  quantizer 
is  a  modification  of  condition  (a) ; 


Re  F (Z) 


1  + 


q 


z-i 

z 


-  \  k'q  |  (Z-l) 

f 


2 


^  0 


for  some  non-negative  q,  where  k  is  the  lower  bound  on 
the  derivative  of  4>(a)  and  its  upper  bound  is  infinite. 
The  local  slope  of  the  quantizing  characteristic 
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curve  has  a  lower  bound  of  zero  and  an  upper  bound  of 

I 

infinite  value,  as  required;  the  quantity  k  is  thus  set 
to  be  zero  in  the  above  condition.  For  this  quantizing 
case,  the  condition  becomes; 


Re  F (Z) 


1  + 


q 


Z-l 

Z 


^0 


for  all  points  on  |Z|  =1  and  some  non-negative  q. 

The  discussion  of  absolute  stability  and  compensation 
that  follows  is  based  on  the  last  condition. 


When  the  value  of  k  is  fixed,  say  two  in  our  case, 

Z-l 


a  minimum  value  must  be  found  for  Re  F(Z) 


1  +  q 


as  Z 


takes  on  values  along  |Z|  =1  and  the  non-negative  q  also 

changes.  An  attempt  to  obtain  analytically  the  optimal 

value  for  q  is  quite  tedious  and  lengthy  as  the  linear  part, 

F(Z),  increases  its  complexity.  A  digital  program  to  find 

the  optimal  q  is  given  in  appendix  (2)  for  a  fairly  general 

form  of  F(Z).  In  this  case  using  the  digital  computer,  the 

optimal  q  has  been  found  to  be  q  =  0.870  which  yields; 

Z-l 


min  Re  F(Z) 


1  +  0.87 


=  -0.7718.  Again,  the  same 
graphical  method  as  developed  in  section  five  and  six  is 
used  here  for  obtaining  stability  and  compensation.  Re¬ 
ferring  to  figure  38,  with  the  aid  of  the  Tsypkin  curve, 
the  minimum  real  part  is  found  to  be  -0.7715.  The  maximum 
gain  allowable  (for  k  =  2)  is  0.648  which  almost  doubles 
the  previous  value  of  1/3  in  Tsypkin' s  case  (q  =  0).  If 
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a  gain  of  K  =  1.16  is  required,  it  is  required  that 

min  Re  F c ( Z )  = (YY~CT~ 16 Y  ~  0.43  2.  To  be  safer,  the 
Tsypkin  curve  is  shifted  down  cutting  the  180°  axis  at 
approximately  -7.5  db  =  0.4265  (180°)  as  shown  in  figure  38. 

In  this  case,  compensation  is  expected  to  be  more 
difficult  due  to  the  fact  that  the  original  gain  without 
compensation  is  already  high.  However,  by  trial-and-error , 
using  the  same  technique  as  before,  the  following  optimal 
values  were  selected; 

<t>  =  30°  ,  co  =  2.05 

m  *  m 


The  compensated  curve,  also  plotted  in  figure  38,  shows 

that  the  resultant  curve  is  below  the  specified  Tsypkin 

curve  for  all  frequencies.  The  minimum  real  part  obtained 

on  the  digital  computer  is  -0.4267. 

For  the  above  values  of  <t>  and  GO  chosen,  locations 

mm  ’ 

of  pole  and  zero  are  given  as;  A  =  -0.4807732,  B  =  0.025311336 
so  that  s 


D(Z) 


=  Z  -  0.025311336 
Z  -  0.4807732 


The  procedure  for  synthesis  is  no  different  from  the  previous 
one.  Since  we  have  a  negative  pole,  the  only  possible  way 
of  realization  is  by  a  digital  computer  as  shown  below  in 
figure  39. 
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Figure  39.  Realization  by  a  digital  computer. 
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8  Conclusion  of  the  Graphical  Technique 

for  Stability  and  Compensation 

The  objective  has  been  to  present  a  new  graphical 
procedure  for  determining  absolute  stability  in  the  large 
and  compensation  of  Tsypkin's  case  and  the  modified  case 
by  Jury  and  Lee.  In  particular,  for  the  compensation  with 
respect  to  the  gain  of  the  linear  plant,  this  method  is 
quite  satisfactory.  No  attempt  was  made  to  obtain  it 
analytically,  due  to  the  mathematical  difficulty  in  this 
kind  of  system.  All  necessary  information  in  terms  of 
frequency  response  is  converted  into  a  gain-phase  diagram 
and  the  whole  procedure  is  carried  out  in  the  Z-domain 
without  any  further  transformation.  The  results  so  ob¬ 
tained  graphically  agree  very  closely  with  those  worked 
out  on  a  digital  computer.  Physical  realization  of  the  lead- 
network  so  developed  is  shown  to  be  possible.  In  Jury  and 
Lee's  case,  some  programming  work  was  performed  to  obtain 
the  optimal  q  which  yielded  the  least  conservative  result. 
This  method  is  general  as  far  as  the  shape  of  the  transfer 
characteristic  of  the  nonlinearity  is  concerned.  However, 
for  finite  local  slopes,  the  last  term  in  Jury  and  Lee's 
criterion  cannot  be  dropped.  In  that  case,  the  digital 
program  must  be  slightly  modified  to  find  the  optimal  q. 

All  other  steps  remain  essentially  the  same. 
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Referring  to  noise  in  the  practical  system  or  in¬ 
accuracy  of  the  quantizer  constructed,  the  stability  con¬ 
dition  has  only  to  be  slightly  modified  to  take  this  into 
consideration.  For  example,  in  our  case,  for  +  10  m.v. 
of  inaccuracy  of  the  electromechanical  relay  used  the 
critical  gain  is  modified  as  follows; 


the  new  sector;  k  -  n  c  v  tt  =  2.04.  The  critical 

s  0 . b  -  U . U± 

gain  is  accordingly  reduced  by  two  percent. 

Although,  a  second-order  system  with  unity  feedback 
has  been  used  to  serve  the  purpose  of  an  illustrating 
example  throughout  the  discussion,  the  usefulness  of 
the  graphical  method  is  not  necessarily  limited  to  such 
a  particular  system. 
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9  Appendix 


1  Scheme  for  starting  the  system  at  the 

sampling  instant 


In  working  with  the  analog  computer  on  sampled-data 
systems,  a  method  is  employed  to  start  the  system  right 
at  the  sampling  instant  with  an  error  of  the  order  of 
microseconds . 

If  the  input  is  a  step-function,  it  is  easy  to  do 
this  by  the  addition  of  a  ZOH. 


lnP^  L  zoH 


tetep  -fu  nclfbv) 


The  actual  input  to  the  system  starts  only  at  the  instant 
of  sampling.  However,  for  other  time-varient  inputs 
such  as  ramps,  or  sinusoidal  inputs,  a  single  ZOH  is  not 
adequate  for  this  purpose.  An  electronic  comparator  is 
used  in  the  new  scheme  as  shown  in  figure  41. 
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When  the  sum  of  inputs  to  the  electronic  comparator  is 
positive,  the  output  marked  by  a  "+"  sign  supplies  a 
positive  voltage  high  enough  to  close  the  switch  in  the 
amplifier  circuit.  Before  the  manual  switch  is  closed, 
the  output  of  the  ZOH  is  zero;  the  input  to  comparator 
is  -5  v.;  the  sign  output  is  negative;  the  switch 

in  the  amplifier  is  open  and  no  input  to  the  system  is 
present.  When  the  manual  switch  is  closed,  the  output 
of  the  ZOH  is  not  +10  v.  until  the  sampling  instant,  at 
which  the  net  input  to  the  comparator  is  +5  v.  When  the 
switch  in  the  amplifier  is  closed,  this  gives  whatever 
input  is  applied  to  the  system.  Once  the  manual  switch 
is  closed,  the  switch  in  the  amplifier  remains  closed. 


*  MOSFET:  metal-oxide  silicon  field-effect  transistor. 
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The  same  principle  can  be  applied  to  control  the  circuit 
generating  the  input  functions,  so  that  everything  can 
be  started  synchronously  at  a  sampling  instant.  Also, 
for  an  integrator,  the  comparator  can  form  a  pair  of 
synchronous  switches  for  reset  and  operate  modes  as 
shown  below  in  figure  42. 


Figure  42.  Electronic  comparator  used  to  control  the 
reset  and  operate  modes. 

Four  electronic  comparators  were  made  together  with 
eight  MOSFET  switches  and  associated  protecting  circuits. 
The  overall  switching  time  is  in  the  order  of  microseconds 
and  the  dead-zone  accuracy  is  +  10  m.v.  No  attempt  was 
made  to  obtain  greater  accuracy  or  speed,  since  the  am¬ 
plifier  noise  on  the  computer  is  about  10  m.v.  and  the 
sampling  time  controlled  by  the  pulses  is  one  millisecond. 


I  serf 
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2  Digital  programs 

Throughout  this  work,  several  digital  programs  were 
written.  They  are  intended  as  tools  only.  These  are 
the  following: 

(1)  Digital  program  for  transient  response. 

Using  equations  derived  from  the  state-variable 
method,  No.  1,  2  of  section  3,  the  transient  response  of 
the  autonomous  system,  subject  to  initial  conditions,  was 
manipulated  by  the  digital  computer.  Extension  of  this 
case  to  one  with  any  input  can  be  easily  obtained  by  slight 
modification  of  the  program.  In  the  actual  program,  using 
computer  variables,  initial  conditions  are  arbitrarily 
chosen  to  be  x  =  3,  y  =  2;  with  a  gain  of  C  =  0.6.  Nine 
intersampling  values  are  calculated.  The  complete  program 
is  shown  as  follows?  • 


, 
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345008  CHfcNG  FORTRAN  SOURCE  LIST- 

ISN  SOURCE  STATEMENT 

0  $IBFTC  S-T  NODECK 

1  DIMENSION  X (200) .Y(200) 

2  C= .  6 

3  T=l. 

4  X(l)=3. 

5  Y( 1) =2 . 

6  K=1 

7  7  H=-X(K) 

10  IF  (H.LT.O.)  GO  TO  4 

13  M=tfK5 

14  GO  TO  5 

15  4  M=H- . 5 

16  5  0=M 

17  300  U=0/H 

20  A= . 1 

21  6  YY=-U*C* (l.-EXP (-A) )*  X(K)+Y(K) /EXP (A) 

22  XX=(1.-U*C* (A-l ,+EXP (-A) ) ) *X(K) + ( 1 . -EXP (-A) )*Y(K) 

23  WRITE  (6,900)  XX, YY 

'24  A»A+.l 

25  IF  (A.GT.l.)  GO  TO  400 

30  GO  TO  6 

31  400  Y(K+1)=-U*C*(1.-EXP(-T) ) *X(K)+Y(K)/EXP (T) 

32  X(K+1)~(1.-U*C*(T-1.+EXP(-T) ) ) *X(K)+(1 i-EXP (-T) ) 
*Y(K) 

33  WRITE  (6,900)  X(K+1) ,  Y(K+1) ,  0 

34  K=K+1 

35  IF  (K.GT.50)  CALL  EXIT 

40  GO  TO  7 

41  900  FORMAT  (1X,4E16.8) 

42  END 


' 
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(2)  Finding  the  optimal  q 

The  following  program  deals  with  the  search 

for  an  optimal  q  for  Jury  and  Lee's  case.  The  value  of 

.  -4  3 

q  is  evaluated  in  the  interval  between  10  to  10  .  The 

linear  system  used  is  quite  general  with  the  form  of; 

G  Zn(Z-A)  (Z-B)  (Z-C)  (Z-D)  (1  +  q 

F  'Z^  ~  {^o)'(z^pr(^)T^sT^u+  j V )  ( Z  - U-  j V ) 

Various  values  of  the  real  constant,  G,  n,  A,  B,  C,  D,  O,  P, 
R,  S,  U,  V  can  be  chosen  at  the  beginning  to  fit  the 
particular  system.  For  example,  in  this  case,  G  =  0.368, 

A  =  “0.717,  B“C=D=0,  0  =  1.0,  P  =  0.368,  R  =  S  =  U  = 

V  =  0,  therefore  n  =  1  was  chosen,  so  that  the  general 
form  becomes; 


F  (Z) 


0.368 (Z+0 .717)  n  ,  Z-l 
(Z-l)  ( Z~0 . 368)  1  q  Z 


The  index  n  becomes  Z  in  the  program.  The  complete  pro¬ 
gram  is  as  follows^ 


,  .  0 
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345008 


FORTRAN  SOURCE  LIST 


ISN 


SOURCE  STATEMENT 


0 

1 

2 

3 

4 

5 

6 
7 

10 

11 

12 

13 

14 

15 

16 
17 
20 
21 
22 

23 

24 

25 

26 
27 

30 

31 

32 

33 

34 

35 

36 

37 

40 

41 

42 

43 

44 

45 

46 

47 

50 

51 

52 

53 

54 


$IBFTC  CHENG  NOLIST , NODECK 

G=EXP(-1.) 

Z=l. 

A=(-l.)* (1.-2. *EXP(-1.) ) /EXP (-1.) 

B= .  0 
C= .  0 
D= .  0 
0=1.0 
P=EXP (-1 . ) 

R= .  0 
S= .  0 
U=.0 
V=.0 

SM=- 1000 . 

BM=1000 . 

Q= .  0 

XX=1.0 

STEP=100. 

1  W=XX*. 017453293 
X=COS (W) 

Y=SIN(W) 

AM=SQRT(X**2+Y**2) 

AMZ=AM**Z 

AMA=SQRT( (X-A)**2+Y**2) 

AMB=SQRT( (X-B) **2+Y**2) 

AMC=SQRT ( ( X-C ) *  *  2+Y*  *  2 ) 

AMD=SQRT ( ( X-D ) *  *  2+Y*  *  2 ) 

AM0=SQRT ( (X-0) **  2+Y**2) 

AMP=SQRT ( (X-P) **2+Y**2) 

AMR=SQRT ( (X-R) **2+Y**2) 

AMS=SQRT( (X-S) **2+Y**2) 

AMU=SQRT( (X-U) **2+(V+Y) **2) 

AMV=AMU 

AMQ=(1.+Q)*SQRT( (X-Q/( 1 .+Q) ) **2+Y**2) /AM 

SIZE= (G*AMZ*AMA*AMB*AMC*AMD*AMQ) / (AMO* AMP* AMR* AMS* AMU* AMV) 
ANZ=Z* 57 . 29 5 7 8* AT  AN2(Y,X) 

ANA=57 . 29578*AT  AN2(Y,X-A) 

ANB=57 . 29578*AT  AN2(Y,X-B) 

ANC=57 . 29578*AT  AN2(Y,X-C) 

ADD=57.29578*AT  AN2(Y,X-D) 

AN0=57.29578*AT  AN2(Y,X-0) 

ANP=57.29578*AT  AN2(Y,X-P) 

ANR=57 . 29578*AT  AN2(Y,X-R) 

ANS=57.29578*AT  AN2(Y,X-S) 

ANU=57.29578*AT  AN2(V+Y,X-U) 
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56 

57 

60 

61 

64 

65 

66 

67 

70 

73 

74 

77 

100 

101 

102 

103 

106 

111 

114 

117 

122 

125 

130 

131 

132 

135 

136 

137 

142 

143 

144 

147 

150 

151 

154 

155 

156 

161 

162 

163 

166 

167 

170 

171 

172 

173 

176 

177 
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ANV-57  *  29578*AT  AN2 (Y-V,X-U) 

ANQ=57 . 29578* (AT  AN2 ( Y, (X-Q/ ( 1 .+Q) ) ) -AT  AN2(Y,X)) 
PHASE=ANZ+ANA+ANB+ANC+ADD+ANQ-ANO-ANP-ANR-ANS-ANU-ANV 
REEL=SIZE*CQS (PHASE* .017453293) 

IF  (BM. LE . REEL)  GO  TO  3 

BM“REEL 

QB=Q 

XXX-XX 

3  XX=XX+1.0 

IF  (XX. GE. 180. 5)  GO  TO  4 
GO  TO  1 

4  IF  (SM.GE.BM)  GO  TO  5 
SM=BM 

QA=Q 

5  Q=Q+STEP 
BM=1000. 

IF  (STEP. LE .0.0005)  GO  TO  50 
IF  (STEP. LE. 0.005)  GO  TO  49 
IF  (STEP. LE. 0.05)  GO  TO  48 
IF  (STEP. LE. 0.5)  GO  TO  47 
IF  (STEP. LE. 5.0)  GO  TO  46 
IF  (STEP. LE. 50.)  GO  TO  45 
IF  (Q.GE.1050.)  GO  TO  11 
XX=1.0 
GO  TO  1 

45  IF  (Q.GE.QQ-l-210. )  GO  TO  12 
XX=1 . 0 

GO  TO  1 

46  IF  (Q.GE.QQ4-21. )  GO  TO  13 
XX=1.0 

GO  TO  1 

47  IF  (Q.GE.QQ+2. 1)  GO  TO  14 
XX=1 . 0 

GO  TO  1 

48  IF  (Q.GE.QG4-.21)  GO  TO  15 
XX=1 . 0 

GO  TO  1 

49  IF  (Q  .  GE  .  QCH"  .021)  GO  TO  16 
XX=1.0 

GO  TO  1 

50  IF  ( Q . GE . QQ+ .0021)  GO  TO  17 
XX=1.0 

GO  TO  1 

11  WRITE  (6,100)  QA , SM 

Q=QA“ 100 . 

QQ=Q 

IF  (Q .GE .0.0)  GO  TO  30 
Q=0.0 
30  XX=1. 
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STEP” 10 „ 
SM=-1000. 

BM— 1000 . 

GO  TO  1 

12  WRITE  (6,100) 
Q=QA-10. 

QQ-Q 

IF  (Q.GE.0.0) 
Q“0 . 0 

31  XX* 1. 

STEP-1 . 

SM— “1000  „ 
BM-1QQQ . 

GO  TO  1 

13  WRITE  (6,100) 
Q-QA-1. 

QQ-Q 

IF  (Q  .GE .0.0) 
Q“0 . 0 

32  XX-1. 

STEP— . 1 
SM— “1000 . 

BM— 1000 . 

GO  TO  1 

14  WRITE  (6,100) 
Q-QA- . 1 
QQ-Q 

IF  (Q.GE.0.0) 
Q-0 . 0 

33  XX- 1 . 

STEP—  e  01 
SM— “1000 „ 

BM— 1000 . 

GO  TO  1 

15  WRITE  (6,100) 
Q— QA” . 01 
QQ-Q 

IF  (Q.GE.0.0) 
Q-0 . 0 

34  XX-1. 

STEP-. 001 
SM— “1000 . 

BM— 1000 . 

GO  TO  1 

16  WRITE  (6,100) 
Q=QA“ .001 
QQ-Q 

IF  (Q.GE.0.0) 
Q-0 . 0 

35  XX- 1. 


QA  ,SM 

GO  TO  31 


QA ,  SM 
GO  TO  32 


QA ,  SM 
GO  TO  33 


QA ,  SM 

GO  TO  34 


QA ,  SM 

GO  TO  35 


2  74 

STEP= . 0001 

275 

SM=- 1000 . 

2  76 

BM-1000 . 

277 

GO  TO  1 

300 

17 

WRITE  (6,100)  SM,QA 

301 

100 

FORMAT  (2E20.8) 

302 

CALL  EXIT 

303 

END 

' 
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(3)  To  plot  the  gain-phase  diagram  and  to 
calculate  the  minimum  real  part. 

After  the  optimal  q  is  found,  the  following  pro¬ 
gram  is  concerned  with  obtaining  information  to  enable 
plotting  of  the  compensated  system  in  the  gain-phase  plot 
as  Z  takes  on  values  on  |Z|  =1  from  coT  =  1°  to  toT  =  180° 
in  steps  of  one  degree.  This  program  is  meant  for  q  =  0.87, 

=  2’05>  ^m  =  30°>  or 


FC(Z) 


0.368 (Z+0 .717) 

(Z-l) (Z-0.368) 


(1  +  0.87  (fEf) 


where  the  expressions  for  finding  B,  A  as  determined  by 

the  choice  of  co  and  <t>  is  included  in  the  program. 

mm 

Notation  q,  B,  A,  u)  ,  4>  become  respectively  in  the 
computer  variables,  Q,  BB,  AA,  WM,  FIM. 

The  complete  program  is  given  as  follows; 
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345008 


FORTRAN  SOURCE  LIST 


ISN 


SOURCE  STATEMENT 


0  $IBFTC  CHENG  NOLIST,  NODECK 

1  A= (1.-2. *EXP (-1 . ) ) /EXP (-1 . ) 

2  B=EXP (-1 . ) 

3  Q= .87 

4  H=1 . 

5  SM=1000. 

6  T=1 . 0 

7  WM=2 .05 

10  FIM=30. *0.017453293 

11  BB= ( 1 ,+SIN (FIM-WM*T) ) / ( SIN (FIM)+COS (WM*T) ) 

12  AA= (SIN (FIM) -COS (WM*T) ) / (-1 .+SIN (FIM-WM*T) ) 

13  WRITE  (6,200)  BB ,  AA 

14  200  FORMAT  (2E20.8) 

15  XX=1. 

16  1  X=COS (XX* .017453293) 

17  Y=SIN(XX*. 017453293) 

20  W=XX* .017453293 

21  AMA=SQRT ( (X+A) **2+Y**2) 

22  AMB=SQRT( (X-l.) **2+Y**2) 

23  AMC=SQRT( (X-B) **  2+Y**2) 

24  ANA=57 . 29578*AT  AN2(Y,X+A) 

25  ANB=57 . 29578*AT  AN2(Y,X-1.) 

26  ANC=57 . 29578*AT  AN2(Y,X-B) 

27  AMQ=SQRT ( ( ( 1 .+Q) *X-Q) **2+ ( ( 1 .+Q) *Y) ** 2) 

30  AMZ=SQRT(X**2+Y**2) 

31  ANQ=57 . 29578*AT  AN2 ( ( 1 .+Q) *Y, ( 1 .+Q) *X-Q) 

32  ANZ=57 . 29578*AT  AN2(Y,X) 

33  SIZ=SQRT( (X-BB) **2+Y**  2) /SQRT ( (X-AA) **2+Y**2) 

34  PD=57. 29578* (AT  AN 2 (Y°X-BB) -AT  AN2 (Y,X-AA) ) 

3  5  AMG= ( H*  B* AMA* AMQ*  SIZ) / ( AMB* AMC* AMZ ) 

36  PHASE=ANA-ANB-ANC+PD+ANQ-ANZ 

37  DB=20 . *ALOG10 ( AMG) 

40  EL=AMG*COS (PHASE* .017453293) 

41  EH“AMG*SIN (PHASE* .017453293) 

42  WRITE  (6,100)  XX,W,E1 , EH, PHASE ,DB,Q 

43  IF  (SM.LT.EL)  GO  TO  7 

46  SM=EL 

47  XXX=XX 

50  7  XX=XX+1 

51  IF  (XX.GT.185.)  GO  TO  3 

54  GO  TO  1 

55  3  WRITE  (6,105)  SM,XXX 

56  100  FORMAT  (7E16.4) 

57  105  FORMAT (/5X, 22HTHE  MINIMUM  REAL  PART ,E12 . 4 , 22H, OCCURS 

AN  ANGLE  0  IF, El 5. 4, 9H  DEGREES.///) 

60  CALL  EXIT 

61  END 
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3  Discussion  of  accuracy  of  the  graphical  method 

When  the  gain-phase  plot  of  the  plant  has  a  very 
steep  slope  near  the  right  end  of  the  graph,  difficulty 
of  matching  the  tangent  point  of  the  two  curves  results 
in  a  problem  of  accuracy.  This  difficulty  can  be  de¬ 
creased  by  redrawing  the  graphs  with  an  expanded  phase 
angle  scale.  It  is  best  to  use  a  digital  computer  in 
order  to  find  the  exact  minimum  real  value  of  the  plant. 
As  far  as  compensation  is  concerned,  a  proper  filter  can 
be  found  by  the  graphical  method  so  that  the  tangent 
point  occurs  at  a  higher  phase  angle.  With  the  aid  of 
a  digital  computer  there  is  no  difficulty  in  determining 
whether  the  requirement  of  the  compensated  system  is 
met  or  not  after  the  compensating  network  has  been 
found  by  the  graphical  method. 
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